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Abstract 

We study the n-level spectral correlation functions of classically 
chaotic quantum systems without time-reversal symmetry. Accord- 
ing to Bohigas, Giannoni and Schmit's universality conjecture, it is 
expected that the correlation functions are in agreement with the pre- 
diction of the Circular Unitary Ensemble (CUE) of random matrices. 
A semiclassical resummation formalism allows us to express the corre- 
lation functions as sums over pseudo-orbits. Using an extended version 
of the diagonal approximation on the pseudo-orbit sums, we derive the 
n-level correlation functions identical to the n x n determinantal cor- 
relation functions of the CUE. 
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1 Introduction 



Quantum systems whose classical limit is chaotic display universal spectral 
statistics. Their spectral correlations depend only on the symmetry class of 
the system, and agree with predictions obtained from averaging over ensem- 
bles of random matricesfH El El Sj. For instance, let us suppose that the 
time-reversal invariance of a chaotic system is broken by applying a mag- 
netic field. Then, in the semiclassical limit, it is conjectured that its spectral 
correlation functions are in agreement with the predictions of the Circular 
Unitary Ensemble (CUE) (or the Gaussian Unitary Ensemble (CUE)) of 
random matrices. 

A way to understand the origins for this universality is provided by semi- 
classics. In the semiclassical theory, the 2-level spectral correlation function 
is expressed as a sum over the pairs of periodic orbits. Berry introduced a 
useful scheme called the diagonal approximation [5]. In this scheme, when 
the time-reversal invariance is broken, only the sum over the pairs of iden- 
tical periodic orbits is taken into account. Then one can derive the smooth 
(non-oscillatory) part of the correlation function in agreement with the CUE. 
Shukla extended this scheme to calculate the ra-level spectral correlation func- 
tions and succeeded in deriving asymptotic forms of the Fourier transforms 
in agreement with the CUE[n]. 

The diagonal approximation brought about a great progress in under- 
standing universality. However, it is able to only partially reproduce the 
random matrix predictions, even if we restrict ourselves to the case of broken 
time-reversal symmetry. It is an approximation which yields only the smooth 
parts of the correlation functions and misses the remaining oscillatory parts. 

Refined schemes to reproduce the full predictions have been developed 
in the study of the Riemann zeta function. It is conjectured that the com- 
plex zeros of the Riemann zeta function are mutually correlated in a similar 
way as the energy levels of chaotic quantum systems without time-reversal 
symmetry[7l |Hl O |10]. Correlation functions of the zeros can be written as 
the multiple sums over prime numbers similar to the periodic orbit sums in 
semiclassics. One is thus able to develop a scheme analogous to the semiclas- 
sical periodic orbit theory. Using additional input from number theory, it is 
then possible to access oscillatory contributions as well. As a result, the full 
CUE correlation functions have been reproduced under certain assumptions 
in several works dH IH [131 • 

In the case of chaotic quantum systems, although the problem is in 
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some respect more involved, the analogous question was addressed in [TBI 
\T7\ \TE\ . Following progress in the method of semiclassical diagrammatic 
expansions fl9 \ [201 [211 [22] . Heusler et al. proposed a way to evaluate semiclas- 
sically both oscillatory and smooth parts of the full 2-level correlation func- 
tion, and obtained results agreeing with the random matrix prediction [17]. 
Keating and Miiller recently gave a justification of Heusler et al.'s argument |18j 
The essential idea is to relate the correlation function to a generating function 
(a ratio of spectral determinants) and then make use of an improved ("re- 
summed" ) semiclassical approximation for the latter, the so-called Riemann- 
Siegel lookalike formula established by Berry and Keating[23l [2^ [25] . For 
systems without time- reversal invar iance, a generalization of the diagonal 
approximation to this setting is sufficient to derive the full 2-level correlation 
function. 

In this paper, we apply a generalization of Keating and Miiller's method 
to calculate the n-level correlation functions of chaotic quantum systems 
without time-reversal symmetry. In §2, we develop a generating function 
formalism for the n-level correlation functions, introduce the Riemann-Siegel 
lookalike formula, and put it into the formalism. In §3, an extended ver- 
sion of the diagonal approximation is formulated and the n-level correlation 
functions are calculated in the semiclassical limit. The resulting formulas 
involve sums over several different contributions generalizing the sum over 
the smooth and oscillatory parts for the 2-level correlation function. In §4, 
we verify that these sums are identical to the n x n determinantal formulas 
of the CUE correlation functions known from random-matrix theory (RMT). 
This is done by establishing their agreement with a representation of the 
random- matrix average obtained by Conrey and Snaith|15] (interestingly, 
this representation had originally been developed to facilitate the compari- 
son to number-theoretic rather than semiclassical results). We thus confirm 
that the known partial results based on the diagonal approximation are ex- 
tended to the forms in agreement with the full random matrix predictions. 
The last section is devoted to a brief discussion on the result. 

2 Generating function 

Let us suppose that H denotes the Hamiltonian of a bounded quantum sys- 
tem which is chaotic in the classical limit. We are interested in the distri- 
bution of the energy levels Ej (the eigenvalues of H). The density of these 
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energy levels 

p{E) = Y,KE ~ E,) (2.1) 
may be separated into the smoothed part 



and the fluctuation around it. Here VL{E) is the volume of the energy shell 
in the classical phase space and / > 1 is the number of degrees of freedom. 

We now want to determine the n-level correlation functions, which are 
defined as 

Rn{ei. ■ ■ ■ , e„) = (fl''^^ ^ ''^) ^^'^^ 

and describe the fluctuation of the energy level distribution around the 
smoothed density p. To obtain a smooth function, we take the average (■) 
over the windows of the center energy E and energy differences e^. 

The idea of the generating function formalism is to represent the level 
densities in (12. 3 p through traces of the resolvent, 

(where E^ = E±iK,, and k, is an infinitesimal positive number), and then ex- 
press these traces in terms of derivatives of the spectral determinant A{E) = 
det{E-H), 

1 d A{E) 



Tr- 



E-H deA{E 

Eq. (l2.3p then turns into 

\j = l [<T,=±1 



(2.5) 



e=0 



271 p J I E + ej + iOjKj — H 



-Zn 



(2.6) 



T]=e 



where 



deide2 ■■■den 

(ei,e2,--- ,e„), = (m, ^2, • ■ • , ^n) (2.7) 
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and the generating function Zn is defined as 



\3=\ l<T,=±l 



We now derive a semiclassical approximation for Z„. Using the Gutzwiller's 
trace formula for chaotic systems [26], we can express the trace of the resolvent 

9{E^) = Tr^^T^ (2.9) 
as a sum over classical periodic orbits a 

g{E^) = -g{E+) - FaTae^'^^'^^^"' ■ (2-10) 

a 

Here Fa is the stability amplitude (including the Maslov phase), Sa is the 
classical action and = dSa/(^E is the period of a. The smoothed part of 
the trace resolvent is written as g{E). It follows from fl2.10p that 

A(^+) oc exp g{E')dE' 

oc exp ^-z7riV(E+) - ^ F.e'^"^^'')/^ j 

Here N{E) is the smoothed part of the cumulative energy-level density: it 
satisfies a relation p = dN/dE with the smoothed part of the energy-level 
density. 

Let us expand the exponential function and write (12.111) as a sum over 
pseudo-orbits A (a pseudo-orbit is a set of component periodic orbits): 

A{E+) oc e-''^^^^^) 5^F^(-1)"-^ e*^?A{i?+)/n_ (2.12) 

A 

Here ua is the number of the component orbits a, Sa is the sum of Sa and 
Fa is the product of Fa. The factor Fa also includes the correction to the 
sign factor (— 1)"^, when identical orbit copies are contained in A. We also 
find from (12. lip that the inverse of the spectral determinant is expanded as 

A{E+)-' oc e*-^^^'') Yl e*''^(^'')/^ (2.13) 



f2.111 
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Similar formulas for A{E~) and A{E~)~^ are obtained by complex conjuga- 
tion. 

However these results do not yet incorporate the unitarity of the quantum- 
mechanical time evolution, and thus the fact that the energy levels are real. 
Berry and Keating argued that the unitarity requirement of the quantum 
dynamics leads to an approximation 

-i7TN{E) 



A{E) oc e" 



(2.14) 



J2 FA{-lT^e'^^^^^/^ + c.c. 

A {Ta<Th/2) 

for a real E. This formula is called the Riemann-Siegel lookalike formula [23t 
[211 after a similar expression in the theory of the Riemann zeta function. 
Here the contributions of "long" pseudo-orbits (pseudo-orbits for which the 
sum Ta of the periods of the component orbits is larger than half the Heisen- 
berg time Th = 27ihp{E)) in (12.121) are replaced by the complex conjugate 
of the contribution from the shorter pseudo-orbits. 

Putting these results into fl2.8p . we obtain an expression 



2TTpi 



E 



exp 



=±1 
=±1 



^ n 



a, 



"^B, JrjSB^{E+rij)/h 



Bj {Tb^<Th/2) 



(2.15) 



with F 



(1) 



Fa and F^"^^ 



F\ (an asterisk means a complex conjugate). 
Here the sums over Aj originate from using (12.131) and its complex conjugate 
in the denominator of (12. Sp . whereas the sums over Bj result from applying 
(I2.14P to the numerator. The sums over Tj make sure that both summands 
in (I2.14P are taken into account. 

Most of the terms in the sum over aj and Tj vanish when they are averaged 
over E, due to the highly oscillatory phase factor. Expanding the exponent 
of the phase factor involving N in (I2.15P as 



exp 



exp 



m ^{ajN{E + e^) - TjN{E + r/^)} 



Tj) + mp{E) ^{ajtj - Tj7]j] 



(2.16) 
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we see that such cancellations can be avoided when 

n 

holds. Hereafter we concentrate on the terms satisfying fl2.17p . 



(2.17) 



3 The n-level correlation functions in the semi- 
classical limit 

To proceed, we introduce four sets /, J, K and L so that 



/ = 


{jV. = 1}, 




J = 


{jWj = -1}, 




K = 


{jk. = 1}, 




L = 


{J\r, = -l}. 


(3.1) 



It follows from flZTTj) that 

\I\ = \K\, (3.2) 

where | | is the number of the elements when is a set. With this notation 
the phase factor involving the actions in (12.151) can be written as 



exp 



exp 



: ^ {ct.Sa, {E + e,) + t,Sb, {E + r/,) } /h 



I jei keJ 



+ J2 Sb, {E + V,) - E Sb, {E + r/,) /h 

jGK k£L 



(3.3) 



In the semiclassical limit, this phase factor also oscillates rapidly for most 
choices of pseudo-orbits, meaning that the corresponding summands in (12.151) 
will be averaged to zero. In order to find the dominant contribution, we need 
to choose the terms with nearly vanishing exponents. For that purpose, 
we simply assume that the component orbits in Aj, j E I and 5^, j E K 
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(contributing with a positive sign in fl3.3p ) are the same as those in A^, A; G J 
and Bj., k & L (contributing with a negative sign), neglecting repetitions. 
We call this scheme the extended diagonal approximation, as it is a natural 
extension of Berry's diagonal approximation [5]. 

Within the diagonal approximation we can now drop the upper limits 
Th/2 for the pseudo-orbits in fl2.15p . This is possible because each periodic 
orbit is now a common component of two pseudo-orbits. Hence its stability 
amplitude is coupled with the complex conjugate to form the absolute square 
in fl2.15p . Weighing orbits with the absolute square of their stability ampli- 
tude is sufficient to ensure convergence even without an upper limit on the 
sum of periods (for energies with an arbitrarily small imaginary part) [IB]. 

In the extended diagonal approximation, each pseudo-orbit Aj (j G /) is a 
union of the disjoint sets AjHAk {k G J) and Aj HBk {k E L). Consequently 
we find a decomposition 

A, 

X 



Analogous decompositions apply to all other pseudo-orbit sums in (12.150 . For 
instance each pseudo-orbit Aj (j G J) is a union of the disjoint sets Aj fl Ak 
{k G /) and Aj f] Bk {k e K). It follows that 

X 



n 



iSA.nB^{E+ej)/h 



J el. 



(3.4) 



n E^. 

kei \AjnAk 



* f^-iSA.nA^(E+ej)/h 
AjDAk^ 



(3.5) 



keK \AjnBk 
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Similarly for j E K and j E L we obtain the decompositions 



fee J \ -BjOAfc 



X 

fceL \Bjr\Bk 



(3.6) 



and 



n J2 ^4n^fe(-l 



^iSjnAj.g-i5'i3jnAj.(-B+r;j)/a 



(3.7) 

If we substitute (13.41) . (13. 5p . (13. 6p and (13. 7p into (12.150 . we see that as antic- 
ipated the actions almost compensate. The only remaining action difference 
is due to the energy arguments being slightly different. This difference can 
be approximated using expansions of the type 

Sa{E + ej) - SAiE + ek) ~ TA{ej - e^). (3.8) 

Our result thus boils down to 



2npi 



X 



n C(-^(e,-r/.)A) n C(-^(r],-e,)A), 
jei,keL jeK,ke,J 

j€i,keJ jeK,keL 
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where the sums over intersections of pseudo-orbits (with the remaining action 
differences fl3.8p ) were written in terms of the dynamical zeta function 



C(.) = ^|F^P(-ire-^^-, 

A 

as)-' = 5^|F^|V^^-. (3.9) 

A 

In order to see the asymptotic behavior of Zn in the semi classical limit, 
we introduce the rescaling 

e, j = l,2,---,n (3.10) 
Znp 

and, noting (12. 2p . utilize the asymptotic formula 

C(s)cxs, s^O, (3.11) 

which holds for chaotic systems [27]. In the semiclassical limit ^ — > 0, we now 
obtain 

Zr, = ( ^ ^ "S^ (-i)\L\Q^{T.,ei^j-j:jej^j-'EjeKVj+j:,eLVj)/2 



3=1 ':>) 



jei j&K 

X ^ . (3.12) 

n^^j^ "^^•) n^^j "^^) 

k&J k&L 

To compare this result to the random-matrix expression in [15], it is 
helpful to adopt a slightly different notation. The sum over all choices for 
the sign factors tj = ±1, j = 1,2, ... ,n is equivalent to summation over all 
ways to write the set {1, 2, ■ ■ ■ , n} as a direct sum of two subsets K and L. 
The direct sum of disjoint sets K and L is defined as the union K U L and 
denoted hy K + L. The corresponding arguments iej, irjj in (13.121) then form 
sets 

A = {iej\jeK}, B = {-iej\j e L}, 

C = {ir]^\jeK}, D = {-zr],\j eL}. (3.13) 
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Moreover the sum over signs aj = ±1 determining the sets / and J can be 
replaced by a sum over subsets 

S = {iej\j e J n K} c A, 

T = {-iej\j ElnL} CB. (3.14) 

Indeed, if we define 

S = A-S, f = B-T, (3.15) 

and = {—a\a E A4} when is a set, the sets of energy increments 
corresponding to / and J can be expressed through S and T as 

S + T- = {ie,\jel}, 

T + S- = {-iejljeJ}. (3.16) 

With these definitions, when we apply the rescahng (13.101) to (12.61) . Eq. (l3.12p 
yields 

K(ei,--- ,en) = Yl (3.17) 

K+L={1,2,- - ,n} 

where 



da 8(3 



■n=e 



and r(y4, B; C, D) is defined as 

r{A,B-C,D) = J2 ^^pUi E E 

SCA,TCB I \Q,g5+T~ /3eT+5- 76C 5eD 

(\SHT\) ^ ^ 

X z{S + T~,f + S';C,D) (3.19) 



with 

z W, X; y, Z) = — . 3.20 

For example, the 2-level correlation function is calculated as 

i?2(ei,e2) = g(0; {-iei, -262}) + g({2ei}; {-^62}) 

+ g({2e2};{-zei}) + g({2ei,2e2};0), (3.21) 
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where 



g(0; {-iei, -zes}) 



9(0ei};{-ie2}) 



9({«e2};{-«ei}) 
g({zei,ze2};0) 



dtidt2 



r(0, {-261, -iie2}; 0, {-2r?i, -im]) 

-i(ei+e2-»?l-'72)/2 



deide2 



■n=e 



1 

4' 



(3.22) 



92 



deidt2 
deide2 



{-^£2}; {-^^2}) 



(ei -62) fa -r/2) 



_j_ g«{e2-ei-)?l+»?2)/2, 

1 1 - 6-^(^1-^2) 

4" (ei-e2)2 ^ 

1 1 - e*(^i-^2) 

4 (ei - 62 
1 

4' 



(e2 - ^2)(^i - ei) 
(62 -ei)fa -r/2) 



T/=e 



2 ' 



(3.23) 
(3.24) 
(3.25) 



Putting and (IX^ into (1^1) . we can readily find a 

compact expression 



i?2(ei,e2) = 1 



sin{(6i-62)/2} 
(ei - e2)/2 



(3.26) 



4 Determinant expressions 

In [12], Conrey and Snaith analyzed the Circular Unitary Ensemble (CUE) 
of random matrices. Based on the Ratios Theorem |28[ 129] (see also [30] [3T]) 
on the characteristic polynomials, they established the formulas 

R^^''''\e,,---,er.)= ^~(^'^) (4-1) 

K+L+M={1,2,- ,n} 
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for the scaled n-eigenparameter correlation functions Here the union 

of the three disjoint sets L and M is {1, 2, ■ ■ ■ , n} and 



q{A;B) 



(4.2) 



T7=e 



with 



f{A,B-C,D)= exp|-5^«-5^/3l2(S + r",f + 5-;C,D). 

SCA^TCB I aeS peT ) 

{\S\ = \T\) 

(4.3) 

The definitions of A, B, C, D, S, T and z(W, X; y, Z) are the same as before 

(dan, dsn and ds;^). 

On the other hand, it is well known that the scaled n-eigenparameter 
correlation functions of the CUE have determinant expressions! 



sin{(ej - efc)/2} 



(4.4) 



J j,k=l,2,--- ,n 



Therefore, in order to verify the same determinant expressions for the scaled 
semiclassical n-level correlation functions (13.171) . it is sufficient to prove that 
f l3Tri) is identical to (lO) . 

To do this, we note that the derivatives ^ (a G A), ^ (/? G -B) in (13.181) 
may act either on the phase factor 



exp <, - 



a 



(4.5) 



or on z{S + , T + S^;C,D). We can thus sum over all ways to split A 
into two disjoint subsets Ai and A2, and then let the derivatives with respect 
to the elements of Ai act on whereas the derivatives with respect to the 
elements of A2 act on z. The corresponding sets of indices in K are denoted 
by Ki and K2. Analogously B is divided into subsets Bi (with derivatives 
acting on cj)) and B2 (with derivatives acting on z), and the corresponding 
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sets of indices in L are denoted by Li and L2. This yields 



E 



n 



ii'+L={l,2,--- .n} S'CA.TCB (|S'| = |T|) \aeAi,f3€Bi 
Ai+A2=A,Bi+B2=B 



d_d_^ 

da 8(3 



Now it is important that the derivatives 



(4.6) 



T7=e 



in (14.61) are nonzero only if 



ri=e 



5 C and T C B2 



(4.7) 



(4.8) 



hold. This is because, for each element iej G S, —iej is included in C 
T + S~ and the corresponding element irjj is included in C so that z{S + 
T~ ,T + S^;C, D) defined in (I3.20p contains a factor —iej + ir/j vanishing for 
r] = e. Nonzero contributions arise only if all such terms are eliminated by 
differentiating z with respect to all elements iej G S. We thus need to have 
S* C ^42. Analogous reasoning leads to T C -B2- 

Now let us consider the phase factor 0. Each derivative of (p with respect 
to the elements a E Ai G S and (3 E Bi <Z T leads to a factor |. If we 
subsequently identify 77 = e, the exponent of (p turns into 



1 
2 



E E 




We thus obtain 
-Rn(ei, 



E ■< 



E 

_R'+L={1,2,- ,n} ScA2CA,rC-B2CB 
(|S| = |T|) 



E^-E"- 

ae5 /3er 

(4.9) 



E 



1 



2l^il+|Si 



X e 



n 



aeA2,/36B2 



9 5 

9a 9/3 



z{S + T~,f + S--C,D) 



ri=e 



(4.10) 
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Eq. fl4.1UI) can be further simplified, if we use the fact that the summands 
in fl4.10p do not depend on how the elements of Ai + Bi are distributed among 
Ai and Bi. In particular, we have 



n 



d_d_ 

da d(3 

d d 



ziS + T-,f + S-;C,D] 



■n=e 



n d^dp<^ nA2 + T~,fnB2 + S-;C2, D2) 



aeA2,l3eB2 



,(4.11) 



r]=e 



where 



C2 = {tVj\jeK2}cC, 

D2 = {-tVj\j ^ L2} C D 



(4.12) 



and all sets on the r.h.s. of (14. lip can be shown to exclude Ai and Bi. 

Let us prove fHTT]) . It follows from (HSj) that S = SnAi + SnA2 = 
Ai + SnA2 andf = f n^i + f n^s = Bi + fnB2. Then, using Ci = C-C2 
and Di = D — D2, we obtain 



n 



[a 



n + ^) 



z{S + T-,f + S-;C,D) 



aeAi+SnA2+T- 
5eDi+D2 



/3eBi+TnB2+5- 

7eCi+C2 



n («+^) n (^+^) 

aeAi+S^nA2+T- 1&C1+C2 
l3eB-i+TnB2+S- 6&D1+D2 

WiW2z{S n + T^, T n 52 + C2, D2), (4.13) 



where 



/3e-Bi 
7ec 



(4.14) 



/3eBi 



7GC1 

6&Di 



&&D2 



7GC2 



15 



and 



n («+'^) n (^+^) 



w, = ''''' 



n n ("+^) 

_oeAi aesnA2+T- 
/S&TnB2+S- l3£Bi 



n n(i+» n u ^^'"'"' 



a — ie.j J V P + i^i 

(4.15) 

We can see that Wi is independent of the elements of A2 U B2 and 



(4.16) 



Ty=e 



for a G A2UB2. (To check fl4.16p . note that only the a and /3 in the last line 
of (14.151) may belong to A2 or B2. If we take derivatives with respect to any 
of these variables, one factor in the product turns into —i , ^'~^\-, or —i /a^~^\o 
and vanishes after setting r] = e.) Therefore nonvanishing contributions arise 
only if the derivatives of (14. 131) with respect to the elements of A2 U B2 act 
on ^(5* n ^2 + T", f n ^2 + S-; C2, D2). Then Wi and W2 can be replaced 
by their special values 

1^1!^=,= 1^2^ = 1 (4.17) 

at 77 = e. 

Thus Eg. (14. lip is proven, and it indeed becomes irrelevant on how the 
elements of Ai + Bi are distributed among the two subsets, or on how the 
corresponding indices in M = Ki + Li are distributed among Ki and Li. We 
can thus stop to discriminate between Ki and Li. This means that we trade 
the sums over K, L and over A2 , B2 in (I4.10p for one over K2 + L2 + M = 
{1, 2, . . . , n}. Since there are 21*^' ways to divide M into Ki and Li we then 
need to multiply with 21^^' which cancels the factor 



1 _ 1 

:l^ll + |Bl 
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The n-level correlation functions can thus be written as 



RnUl 



E 



X 



n 



d d 

da 8(3 



K2+L2+M={1,2,- ,n} SCA2,TCB2 
{\S\ = \T\) 

z{SnA2 + T-,fnB2 + S-; C2, D2] 



■n=e 



(4.18) 



which is identical to the CUE n-eigenparameter correlation functions (14. ip . 
The sets K2, L2, A2, B2 correspond to K,L,A,B in (14. ip . The determinant 
expressions 



Rniei, ■■■ ,en) = det 



sin{(ej -efc)/2} " 
(e,-efc)/2 _ 

for the n-level correlation functions are thus verified 



(4.19) 



j,k=l,2,--- ,n 



5 Discussion 

In this paper the spectral correlation functions of chaotic quantum systems 
with broken time-reversal symmetry were investigated. Using the Riemann- 
Siegel lookalike formula and an extended version of Berry's diagonal ap- 
proximation, we semiclassically evaluated the n-level correlation functions 
identical to the n x n determinantal predictions of random matrix theory. 

Although the semiclassical results are exactly in agreement with the ran- 
dom matrix predictions, we admit that there are higher order nonzero terms 
in the semiclassical diagrammatic expansion, which are neglected in the ex- 
tended diagonal approximation. The exact agreement suggests that those 
terms mutually cancel each other. Such a cancellation was diagrammatically 
verified in [TTl [201 Ell [22] for the 2-level correlation function. It should also 
be verified for the general n-level correlation functions in future works. 

As noted in Introduction, if the quantum system is symmetric under time- 
reversal, it belongs to a different universality class. In this case, the semiclas- 
sical argument becomes more difficult, because the higher order terms in the 
diagrammatic expansion are more involved and give a net contribution [T71 
[20l [211 [22] . The random matrix predictions are derived from the correspond- 
ing Circular Orthogonal Ensemble (COE) (or Gaussian Orthogonal Ensemble 
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(GOE)), and the n-level correlation functions have 2nx2n Pfaffian forms [55]. 
It would be interesting if difficulties are overcome and one is able to see how 
Pfaffian forms semiclassically appear. 
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